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ABSTRACT 

The Hall effect arises in a plasma when electrons are able to drift with the magnetic 
field but ions cannot. In a fully-ionized plasma this occurs for frequencies between 
the ion and electron cyclotron frequencies because of the larger ion inertia. Typically 
this frequency range lies well above the frequencies of interest (such as the dynamical 
frequency of the system under consideration) and can be ignored. In a weakly-ionized 
medium, however, the Hall effect arises through a different mechanism - neutral col- 
lisions preferentially decouple ions from the magnetic field. This typically occurs at 
much lower frequencies and the Hall effect may play an important role in the dynamics 
of weakly-ionised systems such as the Earth's ionosphere and protoplanetary discs. 

To clarify the relationship between these mechanisms we develop an approxi- 
mate single-fluid description of a partially ionized plasma that becomes exact in the 
fully-ionized and weakly-ionized limits. Our treatment includes the effects of ohmic, 
ambipolar, and Hall diffusion. We show that the Hall effect is relevant to the dynamics 
of a partially ionized medium when the dynamical frequency exceeds the ratio of ion 
to bulk mass density times the ion-cyclotron frequency, i.e. the Hall frequency. The 
corresponding length scale is inversely proportional to the ion to bulk mass density 
ratio as well as to the ion-Hall beta parameter. In a weakly ionized medium, the crit- 
ical frequency becomes small enough that Hall MHD is an accurate representation of 
the dynamics. More generally, ohmic and ambipolar diffusion may also be important. 

We show that both ambipolar and Hall diffusion depend upon the fractional ion- 
ization of the medium. However, unlike ambipolar diffusion. Hall diffusion may also be 
important in the high fractional ionization limit. The wave properties of a partially- 
ionized medium are investigated in the ambipolar and Hall limits. We show that in 
the ambipolar regime wave damping is dependent on both fractional ionization and 
ion-neutral collision frequencies. In the Hall regime, since the frequency of a whistler 
wave is inversely proportional to the fractional ionization, and bounded by the ion- 
neutral collision frequency it will play an important role in the Earth's ionosphere, 
solar photosphere and astrophysical discs. 
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1 INTRODUCTION 

In ideal magnetohydrodynamics (MHD), ions and electrons 
are both tied to the magnetic field. When ions are decoupled 
from the field and electrons are not, the magnetic field and 
electrons drift together through the ions and the generalised 
Ohm's law is modified by the Hall electric field, proportional 
to J X B, where J and B are the current density and mag- 
netic field. This modification of ideal MHD is called Hall 
magnetohydrodynamics . 

Hall MHD plays a crucial role in a variety of astro- 
physical, space and laboratory environments, often provid- 
ing the dominant mechanism for plasma drift against the 
magnetic field, from fiux expulsion in neutron star crusts 
(Goldrich & Reisenegger 1992) to angular momentum trans- 



port in weakly ionized protoplanetary discs (Wardle, 1999; 
Balbus and Terquem, 2001; hereafter W99 and BTOl re- 
spectively). The formation of intensive flux tubes in the so- 
lar atmosphere (Khodachenko & Zaistev 2002) , waves in the 
solar wind (Zhelyankov et al 1996; Miteva et al. 2003), prop- 
agation of whistlers in the Earth's ionosphere (Aburjania at 
al. 2005) and sub-Alfvenic plasma expansion (Huba 1995) 
are but a few examples where Hall MHD appears to play 
significant role. In fusion plasmas, the Hall effect can play 
an important role in describing various discharge behaviour 
(Kappraff et al. 1981; Wang & Bhattacharjee 1993). For ex- 
ample, it can significantly enhance the non-Ohmic current 
drive in tokamaks (Pandey et al. 1995). 

Two mechanisms may decouple the ions from the mag- 
netic field under different physical conditions. This has led 
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to distinct approaches being adopted to investigate the role 
of the Hall effect in the dynamics of laboratory (Kappraff et 
al. 1981; Wang & Bhattacharjce 1993; Pandey et al. 1995), 
space (Huba 1995, 2003; Aburjania at al. 2005; Richmond 
& Thayer 2000; Zhclyankov ct al 1996) and astrophysical 
(W99, BTOl, Goldreich and Reisenegger, 1992) plasmas. 

In a highly ionized plasma the Hall effect arises because 
of the difference in electron and ion inertia: ions arc unable 
to follow magnetic fluctuations at frequencies higher than 
their cyclotron frequency, whereas electrons remain coupled 
to the magnetic field. The corresponding physical scale, the 
ion skin depth, is typically much smaller than the scale of 
the system. In this case the Hall effect has typically been 
incorporated by explicitly including the ion-electron drift in 
the induction equation. 

In a partially ionized plasma the Hall effect may instead 
arise because neutral collisions more easily decouple ions 
from the magnetic field than electrons. In this case, the Hall 
scale can become comparable to the size of the system itself. 
Its effects are typically incorporated through a second-rank 
conductivity tensor appearing in a generalized Ohm's law 
(Cowling 1957; Mitchner and Kruger 1973). 

The Hall dynamics of highly ionized and weakly ionized 
plasmas are similar, but occur on very different frequency 
ranges and spatial scales due to the different mechanisms 
responsible for the underlying symmetry breaking in ion and 
electron dynamics. This has led to some confusion in the 
literature, where estimates of the fully-ionised Hall length 
scale have been applied to the ionized component of partially 
ionized media to conclude that the Hall effect is irrelevant in 
circumstances when it is, in fact, crucial (Huba 1995, 2003; 
Bacciotti et al. 1997; Rudakov 2001). 

The purpose of this paper is to clarify the relationship 
between the fully ionized and weakly ionised limits by devel- 
oping a unified single-fiuid framework for the dynamics of 
plasmas of arbitrary ionisation. Our treatment is of necessity 
approximate in the intermediate case, but has the correct 
behaviour in the highly- or weakly-ionised limits and is not 
strongly limited in applicability in the intermediate ionisa- 
tion regime. This allows us to explore the change of scale in 
the Hall effect in moving from fully to partially ionized plas- 
mas and gain a deep physical understanding of the nature 
of the transition between the two ionisation regimes. Fur- 
thermore, this formulation is useful in gaining insight into 
the behaviour of plasmas that are neither fully ionized nor 
weakly ionized (e.g. near a tokamak wall or the surface of 
a white dwarf), when neutral collisions and ionized plasma 
inertia may both be important. 

The paper is organised in the following fashion. In sec- 
tion 2 we derive a set of fluid equation in the bulk frame 
suitable for the weakly ionized medium and the characteris- 
tic scales on which the Hall effect manifests, are discussed. 
In section 3, waves in a partially ionized plasma arc de- 
scribed and the dependence of the wave damping on frac- 
tional ionization in the ambipolar regime is discussed. The 
very low frequency ion-cyclotron and high frequency colli- 
sional whistler is shown to be the two branches in the Hall 
regime. In section 4 wc discuss the potential wide ranging 
applications of this work to laboratory, space, and astro- 
physical plasmas. A brief summary of the results is given in 
the final section. 



2 FORMULATION 

Space and astrophysical plasmas are generally partially ion- 
ized consisting of electrons, ions, neutrals, and charged and 
neutral dust grains. We shall neglect grains in the present 
formulation and consider a partially ionized plasma consist- 
ing of electrons, ions, and neutrals. The dynamics of such a 
plasma is complex but depending upon the physical condi- 
tions pertaining to the problem at hand, reasonable simpli- 
fying assumptions can be made. For example, the dynamics 
of a protoplanetary disc has been investigated by assuming 
that the neutrals provide the inertia of the bulk fluid and 
plasma particles carry the current (W99). This approach is 
reasonable as in a cold protostellar disc, the ionization frac- 
tion (i.e. the ratio of electron to the neutral number density) 
is very low (~ 10~* — 10~^^) and the relative drift between 
ions and neutrals arc small. Therefore, such a description is 
not only economical but also captures the essential physics 
of the protoplanetary discs. However, the inertia of the ion- 
ized components may in general play an important role, e.g. 
near the wall of a tokamak, in the lower part of Earth's F- 
region, at the base of the solar chromosphere, in the outer 
part of AGN discs, in the discs around the dwarf novae etc., 
when the ionization fraction is small and yet not negligi- 
ble. In neutron star crusts too, neutron and proton densities 
are comparable and a multi-component description of the 
strongly magnetized fluid is desirable. In the solar chromo- 
sphere, utilizing three component description, Alfvenwave 
damping have been studied in the context of spicule dy- 
namics (Pontiou & Hacrcndcl 1998). In the solar photo- 
sphere, the effect of ion-neutral damping on the propagation 
of waves has also been recently studied (Kumar & Roberts 
2003). Our aim therefore is to develop an approximate single 
fluid like description of a multi-component, partially ionized 
plasma and demand that it reduces to the fully and weakly- 
ionized descriptions in different fractional ionization limits. 
This approximate formulation will permit us to explore the 
relationship between the onset of the Hall effect due to the 
ion inertia or due to the ion-neutral collisions. Furthermore, 
such a description will provide us the freedom to investigate 
the effect of fractional ionization in various limits on the 
MHD wave modes. 



2.1 A single-fluid model for partially-ionized 
pleisma 

We start with the three-component (ions, electrons and neu- 
trals) description of a partially ionized plasma and reduce it 
to a single fluid description. The continuity equation is 

^+V.{p,vj) = 0, (1) 

where pj — rrij rij is the mass density, Wj is the velocity, and 
Tij and rrij are the number density and particle mass of the 
various components for j = i, e, n. We shall assume that 
the ions are singly charged and adopt charge neutrality, so 
that rii = Ug. The momentum equations for the electrons, 
ions and neutrals are 

j=i,n 
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The electron and ion momentum equations (2)-(3) contain 
on the right hand side pressure gradient, Lorentz force and 
coUisonal momentum exchange terms where Pj is the pres- 
sure, E and B are the electric and magnetic field, c is the 
speed of light, and Vij is the collision frequency for species i 
with species j. The electron-ion collision frequency Pei can 
be expressed in terms of the fractional ionization Xe = ne/rin 
and the plasma temperature Te = Tj = T as 

(5) 



l^ei — 51 Tin T S 



where T and n,i arc in K and cm respectively. The plasma^ 
neutral collision frequency Vjn is 



^jn — Tjn Pn — 



< av >j 
rrin + mj 



(6) 



Here < av >j is the rate coefficient for the momentum trans- 
fer by collision of the particle with the neutrals. The ion- 
neutral and electron-neutral rate coefficients are (Draine et 
al. 1983) 



< av >en 



1.9 X 10"" cm^s"^ 



8.28 X 10" 



'r5 cmS"V 



The density of the bulk fluid is 

P = Pe + Pi + Pn ^ Pi + Pn ■ 

Then deflning the neutral density fraction 



(7) 
(8) 
(9) 



the bulk velocity v = (pi Vi + pn Vn)/p can be written as 

V = {l-D)Vi+DVn. (10) 

Note that we are implicitly neglecting the electron inertia in 
(10), and therefore in the momentum equation (13) below. 

The continuity equation for the bulk fluid is obtained 
by summing up equation (1) for each species: 



|+V.(p.) = 0. 



(11) 



The momentum equation can be derived by adding equa- 
tions (2) - (4) to obtain 

J X B 



dv „ 



Pipn 



■ VdVd 



(12) 



P J c 

where P = Pf, + Pi + Pn is the total pressure, vd = Vi — Vn 
is the ion-neutral drift velocity, and J = UeC {vi — Ve) is 
the current density. 

Defining va ~ B/y/A-K p as the Alfvcn speed in the bulk 
fiuid and Ca = \fyplp as the acoustic speed, we note that if 
PipnV% <C p^i^A + Cs) then we may neglect the vdVd term 
in Eqn (12) and recover the single-fluid momentum equation 

To derive a criterion for this, we estimate by rewriting 
the ion and neutral equations of motion (3) and (4) as 

(Pi l^in + Pe Van) Vd = -pi ^ - V (Pe + Pi) 



, J X B W-e Ven 
~r ~r « 



and 

{pi Vii 



+ Pe Ven) Vd = Pn —J— + V Pn + 



J, 



(14) 



(15) 



dt e 
respectively. Multiplying Eq. (14) by p„ and Eq. (15) by p 

and then adding 



J X B 

{pi Vin + Pe Ven) Vd = D 1- VPn — D VP 



+ 



Pi Pn 



dVD 
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The term in the square bracket can be neglected if 



pi 



(17) 



Then equation (16) can be written as. 



Vd 



Cpi Vi, 
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VP 
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where 



/3, = ^, 



(19) 



is the ratio of the cyclotron frequency of the j^'^ par- 
ticle ujcj = eB/rujC (where e,B ,mj ,c denots electron 
charge, magnetic fleld, mass and speed of light respectively) 
to the sum of the plasma-plasma, and plasma - neutral, 
fjnCoUision frequencies. For electrons Ve ~ v^n + Vei and 
for ions Vi = Vin + Vie- While writing (18), we have used 
pe Ven <Si pi Vin- lu thc Weakly iouizcd limit, when L) — > 1, 
neglecting plasma pressure terms and assuming /3e ^ 1, 
Eq. (18) reduces to the strong coupling approximation, i.e. 
Ui, Ri (J X S)/(cpii/i„)(Shu 1983). 

Equation (18) implies that for gradients with a 
length scale L, and signal speed s, vd ~ pn (1 -|- 
1/ D f3e)/ {piVin L)- The associated dynamical frequency is 



s/ L, so thc requirement pi pn vd ^ P 



that the vdVd term in (12) can be neglected for dynamical 
frequencies satisfying 



Df3e 



^pi pn \\+D (ie 



(20) 



At higher frequencies the single-fluid approximation (13) 
breaks down. Note that this frequency constraint is much 

weaker in the highly-ionized and weakly-ionized limits, for 
which p ~ p„ (D ^ 1). In the appendix we show that 
Eq. (20) is a conservative bound on the dynamical frequency. 
Further, we also show in the appendix that (17) is implied 
by Eq. (20). 

To obtain an equation for the evolution of thc magnetic 
field, we need to derive an expression for the electric field E 
in terms of the fluid properties to insert into Faraday's law 

f = -cvx^;. (21) 

We start with the electron momentum equation (2), which 
in the zero electron inertia limit yields an expression for the 
electric field in the rest frame of the ions: 



+ - xs = 

c 

where 



VPe 



J J xB 

+ - 

a eerie 



■ Vd 



(22) 
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e^Ue 



rUe {fen + V&i) 



(23) 



is the ohmic conductivity and J is given by Ampere's law, 

J=^VXB. (24) 

It is desirable to have an expression for electric field (22) in 
the bulk fluid frame. To obtain this we use Vi = v + Dvd, 
with eq. (18) for vd- Substituting the result into (22) to 
obtain 



dB 

at 



= Vx 



a Pi Vi, 



Pi 



VPn - VPi 



VPe I X B 



(25) 



where wc have neglected the "Biermann's battery" contri- 
bution from the V Pe/eue term in eq. (22) as well as small 
terms of order Df3i/f3e, which is < 10~^. The right hand 
side of this induction equation has convective, ohmic, Hall 
and ambipolar diffusion terms respectively. We note that 
ambipolar term in a partially ionized plasma includes a con- 
tribution from the pressure gradient terms as well from the 
magnetic stresses. 

The relative importance of the various terms in the in- 
duction equation (25) can be easily estimated. The ratio of 
the Hall [H) and the Ohm (O) terms gives H/O ~ , the 
electron Hall parameter. The ratio between ambipolar {A) 
and Hall {H) terms are A/H ^ f3i. In the weak ionization 
(D — > 1) limit, A/H ~ f5i i.e. ion Hall parameter determines 
the relative importance between the Ambipolar and the Hall 
terms. In a highly ionized plasma, D ~ and, the ambipolar 
effect becomes inconsequential. Unlike ambipolar diffusion. 
Hall diffusion does not disappear in the high fractional ion- 
ization limit. 

The ambipolar diffusion terms in (25) arise from D vd X 
B in the X B term in (22) since ViXB = vxB + 
Dvd X B. The terms due to pressure gradients VP X B 
are negligible compared to the inductive term v X B when 



pi Pn 



(26) 



where Cg is some effective sound speed. We note that for 
D fUe ~ 1, Eq. (20) guarantees (26) when va Cs- In the 
opposite limit, when va > Cs, (26) is not implied by (20). 
Our final induction equation without VP X B term becomes: 



dB 

at 



= Vx 



(vxB) J — J X B 



c 

ATTTjA 



(^J X 



X B 



(27) 



where B = B/B, and the Ohmic (r?), ambipolar (r?^) and 
Hall {r]H) diffusivity are 



4 TTO" in Pi Vi, 



,1lH 



cB 



(28) 



Equation (27) is identical to the known expression for a 
weakly-ionised medium (e.g. Kouigl 1989) apart from the 
appearance of the factor D in qA, which suppresses am- 
bipolar diffusion if the ionisation of the plasma is signifi- 
cant. The dependence of ambipolar term on the factor 
was first noted by Cowling (1957). 



To summarize, the single fluid equations have been de- 
rived neglecting electron inertia in the low frequency limit 
given by (20). Pressure gradient terms have been neglected 
in the induction equation (27), which is valid if inequality 
(26) is satisfied. Then equations (11), (13), and (27) along 
with prescriptions for determining P and Ue describe the 
dynamics of a plasma of arbitrary ionization. For example, 
when the plasma is fully ionized, (i.e. 13 — » 0), r = and 
(11), (13), and, (27) reduces to the fully ionized Hall-resistive 
MHD description. In the other extreme limit D ^ I, the 
equations reduce to those describing weakly ionized MHD 
(W99, BTOl). 



2.2 The Hall scale 

In fully ionised plasmas the Hall effect becomes important 
for frequencies in excess of the ion gyrofrequency. In natural 
systems the associated time scales are usually much shorter 
than those of interest and Hall dynamics can be safely ne- 
glected. However, in partially ionised plasmas the Hall effect 
becomes important on longer length and time scales, and in 
weaJily ionised plasmas these may even become comparable 
to the dynamical time scale of the system. 

This behaviour is easily inferred from the fluid equa- 
tions derived in the previous section. If diffusion is unim- 
portant, the characteristic lengthscale of a gradient in the 
fluid associated with frequency ui is L ^ va/co where va is 
the Alfven speed in the total fluid (not just the ionized com- 
ponent). Then comparing the magnitudes of the advective 
and Hall diffusion terms in the induction equation (27), we 
find that the Hall term becomes important for frequencies 
in excess of the Hall frequency 



eB 



m^c p riH 
where the effective ion mass is 
m* = p/rie- 

The corresponding Hall length scale is 

1/2 



Lh = 



Va 



Si = 



va 



(29) 



(30) 



(31) 



where 5i = VAij^d is the ion skin depth with VAi = 
B / ^/TiTpi as the Alfven speed in the ion fluid. 

The Hall effect arises because through an asymmetry in 
the ability of positive and negative charge carriers to drift 
in response to the instantaneous electric field. In the fully- 
ionized limit, for frequencies ujd ^ ^ <^ce electrons are 
able to attain a drift velocity in instantaneous balance be- 
tween electric, magnetic and coUisonal stresses, whereas the 
inertia of the ions prevents them from doing so^. In the 
single-fluid approximation, the ions are tightly coupled to 
the neutrals by collisions so that they are unable to drift 
through them but must carry them along also. Thus they 
pick up the neutral inertia, gaining an effective mass m*, 
and are unable to fully respond to changes with frequencies 



^ The efEects associated with oj > Wce are absent in our esti- 
mate (29) because electron inertia was explicitly neglected in our 
development of the single-fluid equations. 
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in excess of u>h (cf. Pandey & Wardle 2006a, 2006b). Im- 
plicit in this is the requirement that coUisions are able to 
provide the strong coupling between ions and neutrals, as 
noted in eq. (20). 

The condition cj > ujh implies that the Hall term dom- 
inates the inductive term in (27), but does not guarantee 
that it is the dominant diffusion mechanism. As noted ear- 
lier the ratio of the HaJl and Ohmic diffusion terms is ~ /3e, 
whereas the ratio of the ambipolar and Hall diffusion terms 
is ~ D^/3i, so for Hall diffusion to dominate the other mech- 
anisms, we require 



(32) 



Note that in the wcakly-ionizcd limit we recover the stan- 
dard requirement l3i <^ 1 <^ Ps. In the fully-ionized limit 
D^/3i — > and the first inequality is guaranteed. 

2.3 Magnetovorticity 

In the Hall-dominated regime, if the effective ion mass 
m* = p/ue is constant (in space and time) then the con- 
cept of flux freezing can be generalised to the freezing of 
magnetovorticity 



WM = WifS -I- Vxv 



(33) 



into the fluid flow (see e.g., the review by Polygiannakis & 
Moussas 2000). To show this we take the curl of the mo- 
mentum equation (13) and use the identity (V Xv) X v = 
-iVw^ -I- {v ■ V)v to obtain 

^ ' =Vx{v X (yxv)) + Vx 



dt 



pc 



(34) 



where Vp X VP has been neglected. In the absence of mag- 
netic forces, this is the equation for conservation of the vor- 
ticity Vxw. The magnetic term is directly proportional to 
the Hall term in the induction equation 



— = Vx{v X B)+Vx [ 

at \ eue 



and eliminating this term between the two yields 



dt 



Vx {v X um) 



(35) 



(36) 



which shows that the magnetovorticity is frozen into the 
fluid. 

In the limit |Vxv| ^ ujh, this reduces to magnetic 
flux freezing, because Hall diffusion is not significant in the 
induction equation. In the opposite limit, |VXf | ^ u>h, 
ljm reduces to the usual fluid vorticity Vxv because the 
magnetic term is unimportant in the momentum equation. 



3 WAVES IN A PARTIALLY IONIZED 

MEDIUM 

In this section we examine the wave modes supported by a 
partially ionized plasma satisfying eqs (11), (13) and (27). 
We shall assume a homogeneous, uniform background with 
zero flow and investigate the wave properties of the medium 
in various ionization limits. Wc assume that the medium is 
isothermal, i.e. that P = pci with constant sound speed c^. 
The linearized equations for the perturbations 5p, Sv, SB, 
and SJ are 



^ + V.(p^.) = 0, 
d6v 2 r SJ X B 



dSB 

dt 

and 
SJ = 



= Vx 



{6v X B) -5 J — SJ X B 



+ 



c 



(sj X B^ 



X B 



An 



-VxSB. 



(37) 
(38) 

(39) 
(40) 



Note that we do not need an explicit prescription for pi or 
Ue as their perturbations do not appear in the linearized 

induction equation. Assuming that the perturbations have 
the form exp{iujt — ik ■ x), and using (40), equations (37) 
and (38) become 



(jj Sp — pk • Sv = , 

2 f k ■ Sv 
(jjdv = c. 



(41) 

^_ik.B)SB ^iSB.B)k (^2) 



47rp 



Anp 



We define = — fc^c^, dot equation (42) with k, and 
use k ■ SB = to write 

k Sv^ -^!^{B-SB). 

We see that in the incompressible limit, both Sv and SB 
are transverse to the background magnetic field B. Making 
use of equation (43), equation (42) can be written as 



(43) 



ujSv 



i-K p 



{k-B)SB-{^]{SB-B)k 



(44) 



Defining k ■ B = cos 6 and u}\ = v\, and eliminating Sv 
and k ■ Sv from equation (39) we get an equation in terms 
of SB only 

[lu^ — (y'^ + iriA(^) cos^ 6 — irj k'^ ui] SB = 
(SB-B^ (^B-kcosdJ 

-i TjH k^ oj cos ^ ^fe X SB^ , (45) 

where uia = kvA is the Alfven frequency. After some 
straightforward algebra, following dispersion relation can be 
derived from equation (45) 

\io^ — (^ua + i VA i^) cos^ 6 — iri k'' Lu] x 
{ [w^ — (va + ir]AU)) k^ cos^ 9 — irjk^ lj] 



- sin^ 6 



up' 2 

-r7:VA^ir]AU) 



TjH fe* o)^ cos^ 9 = 0. (46) 



In the following sections, we shall investigate the dispersion 
relation, Eq. (46) in various limits. 



3.1 No Hall limit 

In the absence of Hall (i.e. uj <JC ujh), the last term in disper- 
sion relation (46) can be ignored. We note that the modes 
related to the SB \\ B and to SB || fe x B are mixed. For 
example, when the magnetic field perturbation is parallel to 
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the background field, i.e. SB || B, from equation (45), we 
get following dispersion relation 



oj* -i{'nA+ v) - (cs + v\) k^J^ 
+z (A;^ Cg) { {r)A + J?) }t' uj — %}? v'a cos^ ^} = , 



(47) 



which is the second curly bracket in the dispersion relation 
(46). 

In the absence of dissipation (valid for long wavelength 
fluctuations), when r] and r]A — » 0, the roots of the dispersion 
relation are 



1± 



4 k'^ cos^ 
(A;2ci + a;2)2 



1/2^ 



(48) 



The upper and lower sign of Eq. (48) correspond to the fast 
and slow modes of ideal MHD. 

For a cold, coUisional medium, the dispersion relation 
(47) becomes 

- i (7?A + ?7) w A;^ - = 0. (49) 
The real and imaginary part of the root of equation (49) is 

1/2 



Re\jj\ = ±loa 



0.5 



k {r)A + ri) 



VA 



Im\u)\ = 0.5 fc^ {r)A + rf) 



(50) 



It is well known that the waves are damped in the weakly 
ionized coUisional medium (Tanenbaum & Mintzer 1962; 
Kulsrud & Pearce 1969; Ferriere et al. 1988). The damping of 
the waves is not only dependent on the ion-neutral collision 
frequency but also on the ratio of the neutral to the bulk 
mass densities D (Kumar & Roberts 2003). Eq. (50) sug- 
gest that in the absence of Ohmic dissipation, modes with 
the wavelength larger than 



-^cutoff = 



DVA 



(51) 



can propagate in the medium. Clearly, when D = 1, i.e. 
when the medium is weakly ionized, except for a \/2 factor, 
this expression is same as given by Kulsrud & Pcarcc (1969). 
Thus, D = 1 provides the upper bound on the wavelength 
of the damped mode. With the increase in the fractional 
ionization, the cut off wavelength decreases. 

When magnetic perturbation is along fe X -B = n sinS, 
then dotting equation (45) with n sin 6, we get following 
dispersion relation 



(J^ — i (rjA cos^ S + v) 



2 2 

■ lja cos ( 



= 0. 



(52) 



which corresponds to the first square bracket in equation 
(46). The real and imaginary part of the root of equation 
(49) is 



-Re [a;] = iwA 



0.5 



k (rjA cos^ + v) 



VA cos 6 



1/2 



Im[u}] = 0.5 k^ {r]A cos^ ^ + ??) 



(53) 



The normal mode behaviour of the waves are similar in both 
SB ■ B and SB ■ k X B cases except for the cos d reduction 
factor in the later case. We note that equation (53) is iden- 
tical to equation (17) of Desch (2004). The damping of the 



magnetic fluctuations along n is reduced by the cos 8 factor 
in the transverse direction. 

In the limit ^ oo, we obtain the dispersion relation 
found by Desch (2004)) in the Boussinesq approximation 
{dp = 0,SP ^0). 



3.2 Hall limit 

This dispersion relation (46), acquires a familiar form (cf. 
Wardle & Ng 1999, Eq. 25) when wave is propagating along 
the ambient magnetic field {6 = 0) 

up — irfr k^u} — = ^Vh ^ i (54) 

where ?7t = J^a + t]- In the low frequency limit uj <C wa, 
neglecting uJ^ in (54), we get 

Re\iJ\ = ±ujh o , 

Im[oj] = WH . . (55) 



The modified ion-cyclotron mode, Eq. (55) has very low 
threshold of excitation in a weakly ionized medium since 

Re[ij] ~ iajh ~ 0. The ratio of imaginary and real part of 
the frequency only in the presence of ambipolar diffusion 
(?7 = 0) gives 

which is same as the ratio of ambipolar to Hall term in the 
induction Eq. (27). Above expression can also be written 
in terms of Hall and Pedersen conductivities (Wardle & Ng 
1999). We note that when /3i <C 1, damping of the waves will 
be insignificant and system can support very low frequency 
ion-cyclotron modes. 

In the presence of Ohmic diffusion only, above ratio is 



(57) 



/m[w] _ 1 

Re[uj] ~ Te ■ 

Recall that H/O = Pe in the induction equation Eq. (27). 
Therefore when /3e S> 1, i.e. when Hall dominates Ohm, the 
damping of the ion-cyclotron mode is insignificant. We may 
conclude that in the Hall regime, i.e. when 1)2 < K /J^, 
weakly ionized plasma can easily excite very low frequency 
ion-cyclotron mode which will propagate undamped in the 
medium. 

Since the excitation threshold of modified ion-cyclotron 
wave is close to zero this mode will always be present in the 
medium, except when the direction of wave propagation is 
almost transverse to the ambient magnetic field, since for 
oblique propagation iie[a;] ~ ujh cos 6. We note that the 
excitation of the very low frequency modified ion cyclotron 
mode in a weakly ionized medium is a novel feature of the 
Hall MHD. This feature makes it different from a highly 
ionized case. Therefore, in the weakly ionized medium such 
as dark clouds and protoplanetary discs, where luh ~ 0, the 
ion-cyclotron mode is likely to exist in the medium. 

In the high frequency limit uja -C oj, neglecting a;^ in 
(54), and assuming 77 = 0, we get 

IjJh 
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Im[u;] = , 



(58) 



and the ratio of imaginary and real part of the frequency 
gives (3i implying that, the system excites low frequency 
ion-cyclotron and high frequency whistler waves in the sys- 
tem when /3i < 1. 

We note that the nature of the whistler wave in a par- 
tially ionized medium is different from the whistler in a fully 
ionized medium. In ujh —> limit, the whistler frequency can 
become very high but the present single fluid description is 
valid only for whistler frequencies satisfying equation (20). 
In terms of wavelength this constraint becomes 



A>27r I ^ 



1/4 



1+Df3e 



1/2 



r]H cos 6 



1/2 



(59) 



Wo note that the above expression provides a lower bound 
on the wavelength. In a medium such as molecular clouds, 
taking D = 1 and calculating ■qu for a mGauss field with 
Ue ~ .01 cm~^, mi = 30 Trip, rrin = 2.35 rrip and Vni = 
2.8 X 10"^^ s"^ for n„ ~ 10® cm"^, we get A > lO" cm. This 
suggests that single fluid description permits the excitation 
of very small wavelength fluctuations in the cloud. We may 
conclude that the weakly ionized interstellar medium is car 
pablo of exciting high frequency, short {k — » oo) whistlers in 
the medium. 

To summarize, both the modifled ion-cyclotron and 
whistler waves correspond to the short wavelength limit of 
Eq. ( 54 ), i.e. ojh <C wa- In the long wavelength [uja <C ojh) 
limit, we get familiar Alfven wave — oJa- 



4 APPLICATIONS 

In this section wc consider the relevance of the Hall effect 
in fusion, space, and astrophysical plEismgis. We do this by 
adopting typical parameters for the plasmas and examining 
the Hall length and time scales, and the relative magnitudes 
of Hall, ambipolar and Ohmic diffusivities. We also discuss 
its likely implications. 



4.1 Fusion plasmas 

It is well known that in the fusion devices Hall can play an 
important role in discharge behaviour such as the sawtooth 
collapse of the tokamak discharge (Wang & Bhattacharjee 
1993), or, in non-Ohmic current drive schemes (Pandey et 
al. 1995). The Hall effect is also potentially important near 
the partially ionized wall region of a tokamak. Since rin/ni ~ 
(10"^ - 10"") near the waU region (Fiilop et al. 2001), it is 
clear from Eq. (29)-(31) that ujh ~ Wci and Lh w 5i, where 
we have assumed rrin = rrii. Clearly, Hall scaling is similar 
in both fully ionized core and partially ionized wall region 
of the tokamak. Therefore, Hall effect near the wall region 
will be important for coci ^ and scale comparable to the 
ion-inertial scale. 

For typical ion densities ~ lO^** cm^'^ and 10 kG field 
in fusion plasmas, assuming nii = nip we get iOci = 10* s~^. 
The ion Alfven speed is vai = 2.23 x 10** cm, and, the ion skin 
depth is 5i = VAi/oJd ~ 2.23 cm. Adopting ~ 10^ cm as the 
major radius of the tokamak plasma, the Alfven frequency is 
uja = Va ~ lO'^s"^. Therefore, the Hall scale ~ few cm. 



4.2 Ionospheric plasmas 

An important question for magnetosphere-ionosphere cou- 
pling is the interaction between the collisionless magneto- 
spheric plasma and the collisional ionospheric plasma. The 
magnetosphere is well described by the ideal or Hall MHD 
equations whereas the ionosphere is described by the fluid 
equations along with the inertialess plasma determining 
the relationship between the current and the electric field 
through a generalized Ohms law. The transition between 
the two regions is not easily facilitated using these different 
approaches. In particular they mask why Hall operates at 
large scales in the ionosphere, and shrinks to the ion-inertial 
scale in the magnetosphere. The unified set of equations pre- 
sented here treats the ionosphere and magnetosphere in the 
same framework and describes the dynamics of the transi- 
tion region in a consistent fashion. The added bonus of this 
approach is that it explains why the Hall scale shrinks as 
one moves from ionosphere to the magnetosphere. 

To illustrate these points and to gauge the relative im- 
portance of ambipolar. Hall, and Ohmic diffusion in the 
lower ionosphere, we present representative neutral mass 
density, collision frequencies (Akasofu & Chapman 1972; 
Song et al. 2001), and the corresponding Hall-beta param- 
eters, the Hall frequency and Hall length scale in Table 
1. Molecular nitrogen and oxygen are the dominant com- 
ponents of the lower atmosphere, thus we have adopted a 
mean neutral mass m„ — 16 m^. Hall and Ohmic diffusion 
arc dominant in the lower E-layer of the ionosphere. Above 
~ 100 km. Hall diffusion is dominant ujh is very low and 
the corresponding Hall scale is very large (Table 1). With 
increasing height, the density ratio pi/p increases and Hall 
length shrinks becoming of the order of ion-inertial scale 
when pi/ p ^ 1 in the magnetosphere. 

We see from Table 1 that since ratio of the ambipolar 
(A) and Hall (H) terms arc, A/H = rjAhH = -D^ ft ~ 1, 
both these diffusion will operate on an equal footing towards 
the upper E-layer (~ 130 km) and lower F-layer (> 150 km), 
whereas Ohmic diffusion will be unimportant. Observations 
of the partially ionized D and E regions close to the lower 
boundary of the Earth's ionosphere (~ 70 — 140 km), reveal 
the permanent presence of ULF waves. In the E-region of 
the ionosphere, these waves have slow and fast components 
with phase velocities between 1 — lOOms"^ and 2— 20kms~^ 
and frequencies between 10~^ — 10~* Hz and 10~* — 10~® Hz 
respectively, with wavelength > 10"* km and a period of vari- 
ation ranging between few days to tens of days (Zhou et al 
1997; Bauer et al. 1995). Day and night time observations 
gives an order of magnitude difference in the phase velocity. 

The slow ULF waves have been identified as 
Alfven waves, which due to presence of neutrals, converts 
to whistler waves in the E-layer (Aburjania at al. 2005). 
We note that in the lower layer of the ionosphere (D, E 
and lower F -layers, ~ 70 — 140 km), the ionization mass 
fraction could be as low as ~ 10"^^ (Table 1). Hall MHD 
is applicable in this region, and for B = 0.3 G and an ion 
mass mi = mQ+ ~ 10~^^ g, the Hall criterion is cj > u)h = 
10~^ Hz rather than the requirement uj ^ Wd ~ 10^ s"^ that 
would hold if the medium were fully ionized. Thus, waves in 
the 10~^ — 10"** Hz range are most likely whistler waves. 
Identifying the observed wave speed (1 — 2 km) (Aburjania 
at al. 2005) with the whistler, we obtain the characteristic 
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Table 1. Mass Density p, the ratio of ion to neutral mass density Pi/p, ion- neutral, and electron- neutral, 
Ven collision frequencies, ratio of the ambipolar to Hall, tja/vh = D"^ Pi and Hall to Ohm, tih/vo = Pa 
diffusivities along with the Hall frequency u!h and Hall scale length Lh is shown in the table for different 
heights pertaining to Earth's lower ionosphere. A 0.3 G magnetic field has been assumed. 



h(km) 


p(gcm '■'■) 


Pilpn 




(Hz) 


Pen (Hz) 


D^Pi 


Pe 


ujh (Hz) 


Lh (km) 


80 


10-8 


10-12 


7 


• 10^ 


10^ 


10-3 


.07 


10-10 


10^ 


100 


10-9 


10-10 


7 


• 10^ 


10^ 


0.7 


70 


10-8 


iqs 


130 


10-10 


lo-'^ 




102 


3 • 103 


1 


103 


10-6 


10-1 


150 


10-11 


10-4 




30 


8 • 102 


3 


10* 


10-2 


10 



wavelength A ~ 10"^ km at 130 km. Therefore, what is being 
identified as slow Alfvenmode converting to whistler could 
be just low frequency whistler mode without any mode con- 
version. 

In fuUy-ionizcd plasmas the Hall cfltoct considerably 
modifies the classical Kelvin-Helmholtz (KH) and Rayleigh- 
Taylor (RT) instabilities on the ion-inertial scale (Talwar & 
Kalra 1965). Both these instabilities grow faster in the pres- 
ence of Hall effect. We anticipate that similar modifications 
to the KH and RT instabilities will be effective on much 
longer scales and that the Hall effect may facilitate the en- 
ergy cascade from large to small scales in the weakly-ionized 
regions of the ionosphere. 



4.3 Solar atmosphere 

The potential role of the Hall effect in the solar atmosphere 
has escaped the attention of solar community owing to the 
confusion about the Hall scaling in partially ionized plasmas. 
The interaction between the partially ionised solar atmo- 
sphere and fully ionized corona has important consequences 
for the wave heating of the corona. MHD waves can be easily 
excited in the medium by e.g. convcctivc gas motions. A tiny 
fraction of this energy carried by the waves to higher altitude 
would suffice to heat the corona to high temperatures (Priest 
1987). However, the dynamics of the weakly ionized pho- 
tosphere is dominated by collisional effects whereas highly 
ionized corona is described by ideal MHD. It is unclear how 
does one makes a transition between the collisional photo- 
sphere and coUisionless corona. The present formulation not 
only allows us to investigate the dynamics of the collisional 
lower photosphere but also allows tating the proper limit to 
the fully ionized coronal region. As we shall see below. Hall 
may indeed become important in the solar photosphere. To 
show this, we give typical collision frequencies and Hall pa- 
rameters for the solar atmosphere in Table 2 for standard 
solar photosphere and chromosphere models (Vernazza et 
al. 1981; Cox 2000). The proton-hydrogen {H+ - H) elas- 
tic collision cross section is temperature-dependent and at 
0.5 eV is typically 2 ■ 10-"cm2 (Krstic & Schultz 1999). 
The electron-hydrogen (^e~ — H) collision cross-section is 
also temperature-dependent and is 3.5 • 10"!^ crn2 at 0.5 eV 
(Bedersen & Kieffeer 1971; Zecca et al. 1996). Table 2 some- 
what underestimates the collision frequencies as charge ex- 
change between ionized and neutral hydrogen in the solar at- 
mosphere has a large collision cross section ~ 5.6 • IQ-i^ cm2 
(Krstic & Schultz 1999). To calculate the Hall parameters, 
we have assumed a magnetic field B = 100 G. We note here 
that in the solar photosphere electron - ion collisions are 
as significant as plasma - neutral collisions. This could be 



accounted for by multiplying the plasma - neutral collision 
frequencies iu Table 2 by two. However, this does not change 
the overall conclusions of this work. 

Ohmic diffusion dominates in the quiet-Sun photo- 
sphere. However, given the uncertainty about the collision 
frequencies and strength of the magnetic field, Ohmic and 
Hall diffusion may be comparable at the base of the chro- 
mosphere. Then the Hall effect may significantly affect the 
excitation and propagation of waves in the photosphere and 
chromosphere. As is seen from Table 2, near the surface of 
the Sun, the Hall effect will operate over few meters whereas 
near the base of the chromospheres, the typical Hall scale 
may be of the order of few hundred kilometres. Wc note 
that for stronger magnetic fields the Hall scale will be con- 
siderably modified. For example, for a 1 kG field, the Hall 
scale ranges from tens of km near the surface of the Sun to 
hundreds of km as approaching the chromosphere. 

Heating of the solar corona by MHD surface waves is a 
plausible mechanism for explaining the high coronal temper- 
ature (Priest 1987) . Waves in the corona are thought to have 
emerged from the lower photosphere, where they may have 
been excited by footpoint motion of the magnetic field lines. 
Alfven waves propagate to the corona and lose energy by res- 
onance damping, and heat the plasma. However, the lower 
solar photosphere is comparatively cold {T ~ 6000 K) and 
weakly ionized, and ion-neutral collisional dynamics may 
play an important role. We see from the values of f3j in Ta- 
ble 2 that except at h = where Ohm probably competes 
with Hall, Hall is the dominant diffusion mechanism in the 
medium. For a field > lO'^ G (a field probably present in the 
network, internetwork, plague and sunspot Berger & Title 
(2001); Domnguez et al. (2003)) even at /i = Hall will be 
the dominant mechanism. Ambipolar diffusion is unimpor- 
tant in the lower photosphere (?; 1000 km), where l3i <^ 1. At 
higher altitudes, typically between (1000 — 3000) km, am- 
bipolar diffusion will be important, above this the role of 
ambipolar diffusion diminishes as the neutral number den- 
sity plummets rapidly. We infer that Hall effect is present 
over the entire solar atmosphere although with shrinking 
spatial scales with increasing altitude and will modify the 
large scale wave motion in the medium. 

Granulations or convcctivc motions arc believed to gen- 
erate Alfven waves in the photosphere. The Alfven wave is 
a promising candidate for heating and acceleration of the 
solar plasma from coronal holes. High-frequency (~ lO^'Hz) 
ion-cyclotron waves have been proposed as a candidate for 
preferential heating of the heavy ions (Kohl et al. 1998). 
The power spectra of horizontal photospheric motions sug- 
gest that waves with frequencies {10~^ — 0.1 Hz) are present 
at a few solar radii (Cranmer & Ballegooijen 2005). Thus, it 
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Table 2. Same as in Table 1, but for the solar atmosphere. A 100 G magnetic field has been assumed. 



h(km) p(gcm 3) pi/p„ (Hz) i/gn (Hz) ft Pe oj^f (Hz) L^f (km) 

2.77- 10"'^ 10^"' 1.6-10® 1.3 -101° 10""' 10"! 10 10"! 

525 4.87 ■ 10"'^ 10-"' 2.2 ■ lO'' 2 • 10* 10"^ i iq iq 

1000 5.07 • 10"" 10-3 2.2 -10^ 2 • 10*^ 1 30 10^ 10 



should be possible to observe 1 Hz waves. Whistlers may 

be excited close to the footpoint of the flux tube, and 
with decreasing neutral density, this wave will turn into an 
Alfvcnwavc that can propagate to the corona and heat it. 
As the damping of the mode decreases with the increas- 
ing ionization, the waves can propagate almost undamped 
right up to the corona. Thus, in the lower part of the solar 
corona, Hall effect may generate whistler waves. The Hall 
effect operates on scales Lh in excess of a few tens of kms, 
the wave modes due to Hall will be important to various 
coronal heating models. In the fully-ionized solar wind, the 
role of Hall term is constrained by the srnallness of the ion 
skin depth although even then the Hall effect appears to 
significantly modify the surface wave properties in such a 
plasma (Zhelyankov et al 1996). 

4.4 Protoplanetary Discs 

The role of Hall drift on the dynamics of protoplanetary 
discs (PPDs) has been investigated in several papers (W99, 
BTOl, Sano & Stone (2002a,b); Salmeron & Wardle (2003, 
2005); Pandey & Wardle (2006a)). Table 3 gives the collision 
frequencies, plasma Hall parameters and Hall length and 
time scales in a protostellar disc, for a minimum-solar-mass 
nebula model at 5 AU from the central star in the presence 
of a 10-2 G field, using the ionization fraction calculated by 
Wardle (2007). For these parameters. Hall diffusion is im- 
portant for u! > uiH ~ 10~^ — 10-* s-^ comparable to the 
orbital frequency. Near the mid-plane of the PPDs Hall dif- 
fusion will be important whereas towards the surface of the 
disc, ambipolar diffusion becomes dominant. The Hall scale 
Lh ~ 10® km is comparable to tye disc thickness, uggesting 
that Hall operates over the large part of the disc. 



5 SUMMARY 

WeaJily and fully ionized plasmas often exist side by side in 
nature. The transition from weakly to fully ionized plasma is 

typically not abrupt but occurs smoothly, and the transition 
from one to the other has posed a considerable theoretical 
difficulty, requiring separate treatments of Ohm's law for 
weakly (with zero plasma inertia) and fully ionized plasmas. 
In this paper we have developed a consistent MHD frame- 
work for fully, partially, and weakly ionized plasmas. Our 
main results are as follows: 

(i) Our single fluid description encompasses the continu- 
ity, momentum and induction equations (11), (13), and (27), 
along with simple expressions for the Ohmic, Hall and am- 
bipolar diffusivities specified by eqs. (28). These, along with 
prescriptions for determining P and n^, specify the dynam- 
ics of a partially ionized plasma. The equations neglect elec- 
tron inertia limit, and are restricted to low frequencies, i.e. 

^ ^ y/Pn/Pil^el^ni/ (1 + D(5e)- 



(ii) In a partially ionized plasma the frequency ujh above 
which the Hall effect becomes important and the correspond- 
ing spatial scale Lh depend on the fractional ionization, i.e. 



I^H ~ I^ci 

P 




This occurs because the ions are effectively coupled to the 
neutrals by collisions, so that the effective ion mass is 
(p/pi) mi. This behaviour explains why the Hall scale is 
negligibly small in fully ionized plasmas and yet is large 
in protoplanetary discs, and the Earth's ionosphere. It also 
predicts that the Hall effect plays an important role in the 
dynamics of the solar photosphere. 

(iii) We derived a general dispersion relation and showed 
that when ambipolar dominates both Hall [D^ Pi ^ 1) and 
Ohmic {D^ f3e (3i » 1) diffusion, only waves with wavelength 
larger than a cutoff value (Eq. (51)) can propagate in the 
medium. This cutoff wavelength depends on the ratio Pn/ p- 
Low frequency modified ion-cyclotron modes (w = uh) and 
high-frequency whistler [uj — uja/'j.ih) arc the normal modes 
of the medium in the Hall regime. A weakly ionized medium, 
with /3i << 1 will always support these waves owing to 
negligible damping (~ << 1). 

Finally, we emphasise that the single-fiuid formulation here 
will prove useful in studying the coupling between the 
Earth's ionosphere and magnctosphere and also between the 
different solar layers of solar atmosphere. 
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APPENDIX A: 

In this section we derive a general criteria for the validity of 
the single fluid description, i.e. condition under which terms 
^ v'jj in Eq. (12) can bo neglected. Further, wc show that 
the validity condition of the single fluid description ensures 
the validity of equation (17) which allows us to write the 
induction equation (25). We note from Eq. (16) 



Vd 



LOVA 



Dj3e 



(Al) 



and the single fluid description, Eq. (13) is valid provided 
Pi Pn Vo < {v\ + cl) , i.e. 



1 + 



(A2) 

sjpipn V P l^ni J \l + D Pe J ' 

In the following discussion, without loss of generality, we 

shall assume D P^J {1 -\- D P^) ~ 1. The argument is easily 
generalized for arbitrary values of this ratio. It is clear that 
conditiou (A2) is less restrictive than Eq. (20). However, it 
is not cleax if induction Eq. (25) which has been derived by 
assuming an expression for vn, Eq. (18), under condition 
(17) is consistent with Eq. (A2). In order to check the con- 
sistency of condition (A2) with induction Eq. (25), let us 
assume that condition (17) is not valid, i.e. 

UJ . p 



^ni pi 

Writing Eq. (22) in the bulk frame, 

^ ^vuY.B vxB VPe J JxB 

E + D = 1 1 

c c erie cr eerie 



(A3) 



■ Vd ■ 



(A4) 
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We note that if Dvd ^ va then equation (A4) is not valid 
implying invalidity of equation (25). Condition Dvd ^ va 
implies 



> 



(- + -) 

V ^ni Pi J 



(A5) 



Clearly under condition (A5), we axe not allowed to write the 
induction equation (25). We consider two cases: (i) pi > p„ 
and (ii) pi < pn and show that condition (A5) is incompat- 
able with single fluid condition (A2). When pi > pn, it is 
straightforward to see that condition (A5) is never satisfied. 
Thus induction Eq. (25) is valid when under a general single 
fluid condition, Eq. (A2). When pi < Pn, we note that 



and 

U! ^ pi U! 
i^ni Pn ^ni 

Adding Eqs. (A6) and (A7), we get 



l^ni 



l^ni 



(A6) 



(A7) 



(A8) 



Clearly, condition (A3) with pi < p„ implies that Eq. (A2) is 
violated. Thus we have reached a contradiction. Therefore, 
the induction equation (25) is valid under a more general 
condition, Eq. (A2). 

In order to graphically sketch condition (A2), we write 
Eq. (A2) as an equality by multiplying right hand side by 
a small factor /. Assuming / = 0.1 D f3e/ {1 + D f3e), oc = 
\/ pi I pn and y = oj/i/ni, Eq. (A2) becomes 



/Q(l + n -') 
I- fa ■ 

We can rewrite Eq. (A2) in the following form 



/■ 



~ 1_ W V 



(A9) 



(AlO) 



In fig. 1, we plot uj/i^ni against log(pi/p„) using Eq. (A9) for 
/ = 0.1 and / = 0.5. Both in weakly (p, <C Pn), and highly 
(Pn Pi) ionized limits we see from fig. 1 that w ^ i>ni- In 
the weakly ionized limit, ui/uni ~ \/ pn/ pi and thus when 
Pi/ pn Q, uo becomes arbitrary largo (fig. 1). In the highly 
ionized limit, uj/vni ~ 1/ (1 ^ f pi/ Pn)- This results in ar- 
bitrary large w when pn/pi — > /• Therefore, the dynamical 
frequency of a partially ionized plasma is bounded between 
the frequencies of weakly and fully ionized plasmas. 




Figure Al. The numerical solution of tlie equation (A9) for 
D j3e/ (I + D fie) = 1 and 5 corresponding to / = 0.1 and 0.5 re- 
spectively. 



